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It is shown that the Wegner model of disorder contains a system of constraints which is important 
whenever the disorder is not weak and which is responsible for localization in D > 2 dimensions. 
When the disorder is strong the constraints divide the phase space into an exponentially large 
number of domains; the localized phase is a glass. This work is based on recently developed field 
theories which are exactly equivalent to the Wegner model and to an easier variant. 

PACS numbers: 72.15.Rn, 64.60.Cn, 64.60.De, 64.70. Q- 



Disorder, meaning complicated non-repeating struc- 
ture at small scales, is ubiquitous in materials of natural 
origin or made with imperfect human control. Exam- 
ples include sponges, clouds, cellular components, and 
solids containing impurities. Weak disorder means that 
the fine scale structure is close to being a repeating pat- 
tern, so that the material at first examination resembles 
a lattice or a homogenous continuum. Strong disorder 
means that no pattern or repetition can be discerned 
easily. Anderson pointed out that disorder can play a 
crucial role in regulating conduction of electrons or other 
quantum mechanical entities In a disordered mate- 
rial electrons alternate between straight unimpeded mo- 
tion and bouncing off of impurities. If they get turned 
around completely and recross their own path the rules 
of quantum mechanical interference can multiply or can- 
cel the probability of this particular path. As a result 
the electrons may be trapped, unable to move long dis- 
tances, in which case the material does not conduct; this 
is Anderson localization. 

Weakly disordered materials do not conduct in _D = 
{1,2} dimensions, and do conduct otherwise^. Within 
field theory models of disorder, this physics corresponds 
to spontaneous symmetry breaking; the SSB phase 
exhibits long distance correlations describing conduc- 
tion, while the non-SSB phase has no long distance 
correlations [1, Stronger disorder can cause a tran- 
sition from conduction to localization in D > 2 di- 
mensions. The accepted supersymmetric field theory of 
disorder is mathematically unwieldy when the disor- 
der is not weak. This letter analyzes stronger disorders 
using two recently introduced theories which are written 
in terms of two 2x2 matrices, in contrast to SUSY's 
single 4x4 matrix. In effect half of the SUSY theory has 
been integrated, giving a clear mathematical picture of 
intermediate and strong disorders. 

The Models — The standard reference point for stud- 
ies of conduction and localization is Wegner's model of 
weakly coupled grains [1]. Fyodorov showed that the 
zero-dimensional Wegner model is exactly equivalent to 
a field theory with two 2x2 matrices |6|. Disertori ex- 
tended Fyodorov's transformation to any dimension but 
used a simplified variant of the Wegner model[7i]- The 



Wegner model in any dimension is transformed in Ref. 
0. Both Wegner's and Disertori's models describe non- 
interacting electrons moving through a lattice with V 
sites and N states at each site. The lattice geometry is 
completely encoded in a kinetic operator k which has the 
fc > 0, fc|0) = properties required of any Laplacian. In 
Wegner's model the kinetics are constant and the poten- 
tial is random; (H) = efc, {{H — (H))^) cx N^^. e controls 
the balance between kinetics and disorder. In Disertori's 
model both the kinetics and the potential are random; 
{H) = 0, (HH) oc N-^{1 - k). The most important dif- 
ference between the two models is that Wegner's density 
of states p is determined by k, while Disertori's p is the 
semicircular distribution of random matrix theory. As 
usual we study the two point correlator, which can be 
obtained by calculating the second derivative of 



det{E{^ - H)det{El^ - H) 
det{E\^ ~ H)det{E\^ - H) 



and then setting E^ = & . We convert this to a path 
integral using deiA oc J d'tpdipe^'^^'^ and det~^A (x 
/ dSdS*ei^^^^ , and then average over the disorder. 
L = Sign{ ImA). 



Fyodorov converted from the N- vectors S', ■0 to 2 x 2 
matrices[^ [^, 0]. The Hubbard-Stratonovich technique 
converts the vectors to an Hermitian matrix Qf . The 
path integral now depends on dot products Sj ■ Sk, not 
on individual elements of S. Lastly one organizes the two 
S vectors into a, N x 2 matrix and performs the singular 
value decomposition S — WsC. C is 2 x 2 and unitary, 
W is N X 2 and unitary, and s is diagonal and positive. 
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These steps are exact; Disertori's model becomes^, Q 
Z ^ jD [ dQf dQ''L det{Qf - Q''L{1 - k)) 

JQ''>0 

C - -N/2Y,Tr{Q{,Ql + QlLQlL) 

V 

V 

+ {N -2)^Tr\nQlL + Tr\nQl 

V 

- N/2j2{k/{l-k)Uv.Tr{Qf^Qfj 

V1V2 

+ N/2Y,K,..Tr{Ql^LQl^L) (2) 

Wegner's model is exactly equivalent toQ 

Z = jw [ dQf dQ^L dW det(A°(5j,j, - A^) 
Jq''>o 

C = -N/2Y,Tr{Q{,Qf + QlLQiL) 

V 

+ Y,^NTr{QiEf + QlLE'') 

V 

+ [N ~ 2)Tr„, \n{Qf' - lek) - iNeTr^j{SLk) 
+ (A^-2)^rr Ing^L 

V 

= 40 ( A°)~^ ^ r 



-Ell, -B22 are near the band edge or T is quite large. The 
latter condition is avoided if (i?ii — E22)t^P ^ 1 and 
T exhibits SSB. Only in the SSB regime is the saddle 
point spatially uniform. U, T are Goldstone bosons. Dis- 
ertori's saddle point is .Tq — ~ Lp{E) + ts{E), where 
p + is — e"^, sine/) — E/2. When efc <C 1 Wegner's saddle 
point differs by only 0{ek), and when efc ^ 1 the main 
difference is that p is efc's density of states. 

Ward identity — By construction, the Wegner model 
obeys dZ{Ef = E^)/de = 0. The same Ward identity 
governs Disertori's model if we define e = |fc|. We inven- 
tory factors of e, assuming efc ^ 1 for Wegner's but not 
Disertori's model. The W integral's logarithm is dom- 
inated by the self-energy S(s) cx f°N > Nek{s). The 
x-f,x^ integrals. Van der Monde determinants, and re- 
mainder of the Lagrangian are controlled by a;Q,XQ oc 1. 
Therefore any e's produced by the t/, T integrals must be 
cancelled by the Q-^ — determinant. In the SSB phase 
U,T are controlled by the kinetic terms [s'], and produce 
g-2(v-i)^ In the localized phase ?7, T's low-momentum 
behavior should be independent of the kinetics; the inte- 
grals should scale as e~2(v-i-cv)^ where (^V is the vol- 
ume of U, T's phase space which exhibits large fluctua- 
tions. This links the determinant's behavior to localiza- 
tion; it must scale as e^*^^~^~^^-'; the question is how? 

Disertori's determinant is det(aD = — Q'' + Q''k). 
In the SSB phase Wegner's determinant is det(aw w 
~ — lek), plus corrections in powers of U,T,W^s 
fluctuations. Outside of SSB one may expand aw in 
powers of efc <^ 1; ignoring any resonances of Q-^ — lek 
one obtains aw ~ — + 0{ek). This form can be 
rewritten as det(a;'^ — x'' — leky 
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Sy^^^ = Cl^Sy,WlWy,Sy,Cy,, ClslCy = Ql (3) either cancel or add to x"; consequently x^ - x" has two 



?; is a spatial index specifying the lattice site, while 
i,j index the 2x2 Hermitian matrices Q^^Q^. The 
argument of the — Q'' determinant (on the first 
line of both equations) lives in a 2 x 2 x ^ space in- 
dexed by ijv. We consider the Retarded-Advanced 
part of the two point correlator; L = 0-3. We decom- 
pose Qf = Ux^W and Q''L = Tx'^T-^; dQf dQ^L = 
dUdxf dTdx^ n„^ydM(4) ^Idmi^l)- U is unitary. 
T is pseudo-unitary and unbounded, in order to preserve 
the sign signature of Q^Us eigenvalues [5|]. 

Interactions between U and T are mediated only by 
the — determinant. Leaving it aside, the U sector 
is the classical Heisenberg model, while the T sector is an 
hyperbolic sigma model. Both models have been proven 
rigorously to exhibit SSB in _D > 2 dimensions (ol. [lo|. 
The T sector has no phase transition and always exhibits 
ITU . Therefore the determinant is responsible for 
all localization occuring in D > 2 dimensions [lH . 

Saddle point — A separate paper carefully applies the 
saddle point approximation to Disertori's model, and 
partially analyzes the Wegner model[§]. In summary, 
deviations from the saddle point are 0{N~^/'^) unless 



sectors described by two projection operators Pi + Pg 
1. Its large Pi eigenvalues are ±2p, while its small Pg 
eigenvalues are determined by 0(A^^^/^) fluctuations in 
x^ , x^ and i^n— £^22 ^ 1 splitting. If there are no angular 
fluctuations {U — T — 1) a has 2(V^— 1) small eigenvalues 
controlled by ek, 2 smaller ones caused by fcjO) = 0, and 
2F large eigenvalues controlled by ±2/5-|-eA; [8|]. However 
because U, T do fluctuate the kUV'^TT'^ term can mix 
the Pi.Pg sectors, lifting QV eigenvalues to 0{2p) to fulfill 
the Ward identity. 

Optimization Principle — The — determinant is 
the product of 4:V eigenvalues, each of which multiplies 
the entire path integral. The Van der Monde determi- 
nants also multiply by 4V factors. Therefore the path 
integral is governed by an optimization principle: Q-^ 
and always optimize themselves to maximize the de- 
terminants. This principle is limited by the Lagrangian: 
its kinetic terms favor SSB where U, T's fluctuations are 
0{{Nek)~^/'^), while x^ ,x^ fluctuations are restricted to 
0(JV~^/^) plus a response to i7, tQ. In Wegner's model 
when the disorder is small aw is dominated by ek, so 
AV — 2 eigenvalues are pre-optimized without any need 
for mixing, and the kinetics ensure SSB. The remain- 



3 



ing 2 very small eigenvalues caused by fc's zero mode 
fc|0} = cause 0(A^^^/^) changes in x^,x'', and are 
the source of the oscillatory component of Wigner-Dyson 
level repulsion 8]. As the disorder is increased ek's small- 
est eigenvalues become smaller than 2/5 oc (efc)~^; they 
are no longer pre-optimized. The optimization princi- 
ple starts to reward — Pi mixing at small momenta 
ek{s) < 2p and large fluctuations of U,T at large wave- 
lengths. This trend continues into the ek <^ 1 regime, 
where the determinant favors large fluctuations at all 
wavelengths and the kinetic penalty on fluctuations is 
considerably weaker. 

Constraint Structure — Each eigenvalue and Van der 
Monde factor A creates a forbidden hypersurface in 
Qf,Q'"s phase space defined by A = 0. Configurations on 
those hypersurfaces are disallowed; each hypersurface is 
a constraint on Qf,Q^. A logarithmic potential In A pe- 
nalises non-optimized configurations, ensuring that they 
do not contribute to the path integral. In the SSB phase 
only two constraint hypersurfaces lie inside the accessi- 
ble phase space, but as efc's spectrum descends below 2p 
additional constraints are activated. All of the Q-^ — 
constraints intersect the accessible phase space when SSB 
is frustrated at all wavelengths. Leaving aside details of 
global topology and counting, I expect the number of dis- 
joint domains and local minima to be exponentially large, 
at least e*^ but possibly e^^ if the Van der Monde con- 
straints are active. 2{V — 1) of the determinant's eigen- 
values are 0{k <^ 1) when ,Q'' are spatially uniform, 
which implies a near intersection of the constraints and 
prohibits the existence of any spatially uniform local min- 
imum. Fluctuations around each local minimum might 
be small and regulated by the kinetic terms, but the min- 
ima themselves must manifest large fluctuations. 

The constraints' physical meaning is connected with 
the fractional form seen in equation[T] That form, includ- 
ing both the numerator and the denominator, was needed 
to properly average over all disorders, giving each possi- 
ble disorder realization an equal weight. If the disorder 
is small then little reweighting is necessary, but when the 
disorder is large reweighting becomes very important, as 
do the — Q'' constraints. 

Resonances — As fc's lowest eigenvalues sink below 2p, 

the eigenvalues manifest propagating resonances 

which may stimulate U, T to leave the SSB phase. The 

same resonances are manifested also in the sigma model 

controlling [/, T. In the Wegner model these resonances 

are manifested in a;-^,a;'"s Hessian, —NSy^y^ ~ N{ix + 

^k)viV2 (*-^ + ^^)v2Vi ■ During the crossover between small 

and large disorder a;o oc mm(l, {ek)~^) passes through 

the spectrum of efc; the Hessian's kinetic component ex- 

f 

hibits one or more resonances as Xq passes by individual 
eigenvalues of ek. Close to the weakly disordered regime 
the Hessian's kinetics are dominated by small momenta 
and the resonance width is small, but progressively all 
wavelengths become involved. 



Phases and transitions — We imagine starting with 
weak disorder and and later increasing it. Ref. [8| ana- 
lyzes the regime characterized by eko = Eth ^ 1 in the 
Wegner model and by 1 ^ fco ^ N~^/^ in Disertori's 
model, where ko is k's smallest non-zero eigenvalue. In 
this regime mathematical control over Disertori's model 
has been obtained, and a partial analysis suggests that 
the Wegner model might also be controlled. The results 
reproduce and extend standard SUSY results concerning 
the two point correlator, anomalously localized states, 
and localization in D = {1, 2} dimensions. The — 
determinant is nearly static, but its vestiges are responsi- 
ble for Wegner-Dyson statistics. The dominant physics is 
scattering, and is encoded in the sigma model Lagrangian 
controlling C/, T . The Mermin- Wagner theorem I4I says 
that scattering ensures localization in large volumes in 
D — {1, 2} dimensions. The scattering by itself is unable 
to achieve localization in Z) > 2 dimensions, but is able 
to produce individual anomalously localized states. 

As the disorder is increased 2/5 approaches efco and 
xf,x^ develop critical fluctuations at long wavelengths. 
Crossing this threshold activates the first of many new 
constraints, each of which pushes U, T to develop large 
long- wavelength fluctuations with the help of x^ , x'"s res- 
onance. 

A straightforward analysis suggests that the new con- 
straints become active in sequence starting with long 
wavelengths, and that the mobility edge is a function 
of the wavelength, with conduction at small scales and 
localization at large scales. However this is not the last 
word: the model might be able to maximize a few long- 
wavelength constraints without making a full transition 
from 0{{Nek)^^^'^) to 0(1) fluctuations; in this case one 
might see subdiffusive conduction below some threshold 
number of constraints and localization above that thresh- 
old. Alternatively, some collective dynamics might make 
the model unresponsive to new constraints until a thresh- 
old is reached; constraint optimization could occur en 
masse in a single Anderson transition. What is certain 
is that at large enough disorder SSB is completely frus- 
trated with large fluctuations at all wavelengths and all 
of the — Q'' determinant's constraints are active. 

Renormalization group — fc's smallest eigenvalues scale 

as 

therefore increasing the length scale is roughly 
equivalent to increasing the disorder. The renormaliza- 
tion group starts in the weak disorder phase and sees res- 
onances and new constraints as the volume is increased. 
However fc's large eigenvalues are independent of V , so 
the V 00 limit brings only a fraction of fc's eigenvalues 
below 2/5. Anderson localization vs. diffusive or subdif- 
fusive conduction is determined by whether that fraction 
is able to frustrate SSB at long wavelengths. 

We compare this picture to the SUSY picture of 
renormalization [l3| . In the SUSY approach the renor- 
malization group equation's initial value is the SUSY 
sigma model, which is a weak-disorder approximation of 
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the Wegner model. The standard SUSY sigma model re- 
sults (without renormalization) can be reproduced by ap- 
proximating the Q-f — determinant to the point where 
only two constraints are active Q. It would be remark- 
able if the SUSY renormalization group could recover the 
physics of the full Q-^ — determinant from this very 
simplified starting point. 

The renormalization group may be at least partly su- 
perfluous, since the models considered here are exactly 
equivalent to the original Wegner and Disertori mod- 
els including all band and geometric information. They 
avoid the Grassman degrees of freedom found in the 
SUSY model, but by exactness retain their physics; in 
this sense the — determinant should be duplicated, 
not improved upon, by an exact integration of the SUSY 
renormalization group. 

Disertori 's Model — This model gives some hints of an 
Anderson transition around the fco ~ N^^^^ threshold. 
At this point the spectral gap beween k and x-^ — .t^'s 
small eigenvalues disappears, which facilitates mixing, al- 
though the mixing's effects may be limited by N^^^^. 
More interestingly, at the threshold , x^ undergo a 
phase transition as they begin to optimize the Q-^ — 
determinant's eigenvalues. The order parameter is an 
0{N~^/'^) splitting between x^ and x^ enforced by the 
constraints' logarithmic repulsion. A naive computation 
of the eigenvalues' Hessian using this splitting reveals 
a pole on the energy shell k{s) « N~^/'^; x^,x^ orga- 
nize themselves across long distances in order to set an 
0{N^^/^) lower bound on the determinant's eigenval- 
ues. The resonance's length scale becomes progressively 
smaller as k decreases. Its width is 0{N~-^/^), since its 
interaction with U, T is controlled by derivatives of the 
determinant's logarithm, which are proportional to in- 
verse powers of the splitting. 

As mentioned earlier there are exponentially many lo- 
cal minima corresponding to opposite signs of the split- 
ting; determining the resonance's effects on U, T requires 
first integrating x-^,x^. Here lies a mystery. If SSB con- 
trols T's unboundedness then the saddle point integral 
is Gaussian. It is a simple matter to establish that the 
averaged determinant is bounded below by N^'^^ , not 
by the N'^ factor one would expect from the eigenvalue 
splitting; on average the splitting is invisible! Adding to 



the mystery, this disappearing act occurs only in calcu- 
lations of two-point correlators and not in higher order 
correlators, i.e. if one starts with three or more factors in 
equation [T] then the averaged determinant clearly mani- 
fests the eigenvalue splitting. 

If the Anderson transition does not occur at fco « 
7V~^/^, it must occur at fc « iV~^, because after this 
point the kinetic term Nk is unable to regulate fluctua- 
tions. If one pushes further to the point where fc's spec- 
trum is partly or completely negative, the kinetic term 
rewards U, T's fluctations rather than penalizing them. 
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